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We present an empirical analysis of the network formed by the trade relationships between all 
world countries, or World Trade Web (WTW). Each (directed) link is weighted by the amount of 
wealth flowing between two countries, and each country is characterized by the value of its Gross 
Domestic Product (GDP). By analysing a set of year-by-year data covering the time interval 1950- 
2000, we show that the dynamics of all GDP values and the evolution of the WTW (trade flow 
and topology) are tightly coupled. The probability that two countries are connected depends on 
their GDP values, supporting recent theoretical models relating network topology to the presence 
of a 'hidden' variable (or fitness). On the other hand, the topology is shown to determine the GDP 
values due to the exchange between countries. This leads us to a new framework where the fitness 
value is a dynamical variable determining, and at the same time depending on, network topology in 
a continuous feedback. 



I. INTRODUCTION 



The globalization process of the economy is highlight- 
ing the relevance of international interactions between 
world countries. The world economy is evolving towards 
an interconnected system of trading countries which are 
highly heterogeneous in terms of their internal activity. 
As many other socioeconomic and financial systems with 
heterogeneous interacting units, the global economy ex- 
hibits complex structure and dynamics, and can there- 
fore be studied exploiting the tools of modern statistical 
mechanics. 

A role of primary importance is played by the network 
of import/export relationships between world countries, 
or World Trade Web (WTW in the following). Recent 
empirical studies [l|, H, [3, H have focused at the WTW 
as a complex network 0, [y, 0] and investigated its ar- 
chitecture. A range of nontrivial topological properties 
have been detected [3, 0, H, 0] and found to be tightly 
related to the Gross Domestic Product (GDP in the fol- 
lowing) of world countries [H, H[. On the other hand, 
the economic literature has recently dealt with the study 
of the GDP per capita across countries, looking for pat- 
terns and trends displayed by it. Some of these results 
have pointed out that the GDP per capita displays com- 
plex dynamical behaviour and a power-law distribution 
across countries [1, 0] . 

In the present work we extend these analyses and ad- 
dress the above points simultaneously. In particular, we 
are interested in determining empirically the effects that 
the WTW topology and the GDP dynamics have on each 
other. 



II. DATA SET AND DEFINITIONS 

The results of the present work are based on the em- 
pirical analysis of a large data set reporting the an- 
nual values of the population size pi (t) and of the GDP 
per capita Zi(t) of each world country i for each year t 
from 1948 to 2000, together with the annual amount of 
money fij(t) flowing from each country i to each coun- 
try j due to exports from j to i between 1950 and 2000. 
Trade, population and GDP figures are evaluated at the 
end of each year. Since in the following we are inter- 
ested in a simultaneous study of GDP and trade data for 
each country, we restrict our analysis to the time interval 
1950-2000 when both sources of information are avail- 
able. We therefore set our initial time to to — 1950. Our 
main interest is the total economic activity of each coun- 
try, therefore we multiply each GDP per capita Zi(t) by 
the population size pi (t) of the corresponding country to 
obtain the total GDP Wi (t) = Zi (t)pi (t) of that country. 

Trade data are expressed in current U.S. dollars (that 
is, in terms of the value $ t of one U.S. dollar in the re- 
ported year t), while GDP data are available in current as 
well as in 1996 U.S. dollars ($1996)- The use of a standard 
reference money unit such as $1996 factors out the effects 
of inflation and allows a comparison between figures cor- 
responding to different years. Therefore in the following 
we rescale the trade data corresponding to each year t 
to 1996 dollars. This means that Wi(t) will represent the 
real GDP of country i. A curve of the time dependence 
of the ratio r(t) = $t/$igg6 of the value of current U.S. 
dollars to their 1996 reference value is shown in fig. [1] 
In what follows, both fij(t) and Wi(t) will be expressed 
in millions of 1996 U.S. dollars (M$ 1996 ). 

The number N(t) of world countries is monotonically 
increasing in time during the considered time interval, 
and it grows from A(1950) = 86 to TV(2000) = 190. This 
means that the WTW is a growing network. In fig. [5] 
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FIG. 1: Temporal evolution of the ratio r(t) = $t/$i996 from 
t = 1950 to t = 2000. 



activity is therefore given by 

Fi(t) = ft (t) - fr\t) (3) 

We finally dehne the adjacency matrix elements as 



" 1 if /y(t) = 



(4) 



and the in-degree kl"(t) and out-degree k° ut (t) of a coun- 
try i at time t as 



N(t) 

3=1 
iv(t) 

3=1 



(5) 



(6) 



representing the number of countries to which i exports 
and from which i imports respectively. 
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FIG. 2: Temporal evolution of the number of world countries 
N(t) from t = 1950 to t = 2000. 



we plot the time evolution of N{t). The reason for the 
increase of N{t) between 1960 and 1990 is mainly due 
to the declaration of independence of several countries, 
while the steep increase around 1990-1992 is due to the 
origin of many separate states from the former Soviet 
Union. 

The total trade value of exports and imports by i 
to/from all other countries will be denoted by f m {t) and 
f out (t) respectively, and it can be expressed as 



JV(t) 

/:";/; »< j/i 

3=1 
Nit) 

/r*(*) = E/«(*) 

3=1 



(1) 



(2) 



III. GDP: DEFINITION AND EMPIRICAL 
PROPERTIES 

The Gross Domestic Product wi of a country i is de- 
fined as (see for example www.investorwords.com) the 
"total market value of all final goods and services pro- 
duced in a country in a given period, egual to total con- 
sumer, investment and government spending, plus the 
value of exports, minus the value of imports". In other 
words, there are two main terms contributing to the ob- 
served value of the GDP u>j of a country i: an endogenous 
term Li (also konwn as internal demand) determined by 
the internal spending due to the country's economic pro- 
cess and an exogenous term Fi determined by the trade 
flow with other countries. The above definition can then 
be rephrased as 



Wi(t) =Ii{t)+Fi(t) 



(7) 



where Fi(t) is defined by eqs.(T]), © and ©. The above 
definition anticipates that the GDP is strongly affected 
by the structure of the WTW. The characterization of 
the interplay between the GDP dynamics and the WTW 
topology is the main subject of the present work. Be- 
fore addressing this issue in detail, we first report in this 
section some empirical properties of the GDP. 

A fundamental macroeconomic question is: how is the 
GDP distributed across world countries? To address this 
point we consider the distribution of the rescaled quantity 



Xi(t) 



(w(t)) 



(8) 



where (w) = wt/N is the average GDP and u?<r(i) = 

SfcLi w i(t) i s tne total one. In fig. ([3]) we report the 
cumulative distribution 
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(9) 
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FIG. 3: Normalized cumulative distribution of the relative 
GDP Xi(t) = Wi(t) I {w)(t) for all world countries at four dif- 
ferent snapshots. Inset: the same data plotted in terms of 
the rescaled quantity yi(t) = Wi(t)/wT(t) = Xi(t)/N(t) for 
the transition region to a power-law curve. 



for four different years in the time interval considered. 
The right tail of the distribution roughly follows a 
straight line in log- log axes, corresponding to a power- 
law curve 



p>(x) oc x 1 " 



(10) 



with exponent 1 — a = — 1, which indicates a tail in the 
GDP probability distribution p(x) oc x~ a with a = 2. 
This behaviour is qualitatively similar to the power-law 
character of the per capita GDP distribution [H, 0] . 

Moreover, it can be seen that the cumulative distribu- 
tion departs from the power-law behaviour in the small 
x region, and that the value x* where this happens is 
larger as time increases. However, if Xi(t) is rescaled to 



Wj{t) _ Xj(t) 

w T (t) ~ N(t) 



(11) 



then the point y* — x* /N rs 0.003 where the power-law 
tail of the distribution starts is approximately constant 
in time (see inset of fig. [3]) . This suggests that the tem- 
poral change of x* is due to the variation of N(t) affect- 
ing (w(t)) and not to other factors. This is because the 
temporal variation of N(t) affects the average-dependent 
quantities such as x: note that, while for a system with a 
fixed number N of units wt would be simply proportional 
to the average value (w), for our system with a varying 
number of countries the two quantities can be very dif- 
ferent. In particular, the average values of the quantities 
of interest may display sudden jumps due to the steep 
increase of N(i) rather than to genuine variations of the 
quantities themselves. 



IV. EFFECTS OF THE GDP ON THE WTW 

In a recent work [3j it was shown that the topology of 
the WTW, which is encapsulated in its adjacency ma- 
trix a,ij defined in eq.(|l]), strongly depends on the GDP 
values Wi. Indeed, the problem can be mapped onto the 
so-called fitness model [HI, [13] where it is assumed that 
the probability Pij for a link from i to j is a function 
p(xi,Xj) of the values of a fitness variable x assigned to 
each vertex and drawn from a given distribution. The 
importance of this model relies in the possibility to write 
all the expected topological properties of the network 
(whose specification requires in principle the knowledge 
of the N 2 entries of its adjacency matrix) in terms of only 
TV fitness values. Several topological properties includ- 
ing the degree distribution, the degree correlations and 
the clustering hierarchy were shown to be determined by 
the GDP distribution [3J]. Moreover, an additional un- 
derstanding of the WTW as a directed network comes 
from the study of its reciprocity [HI], which represents 
the strong tendency of the network to form pairs of mu- 
tual links pointing in opposite directions between two 
vertices. In this case too, the observed reciprocity struc- 
ture can be traced back to the GDP values [14]. All 
these results were also shown to be robust in time and 
to be displayed by all snapshots of the WTW [H • In this 
section we summarize and further extend these analyses. 

Combining the results presented in refs.[3,[I3.[l3|. the 
probability that at time t a link exists from i to j (a,j = 
1) is empirically found to be 



1 + /3(t)Xi(t)Xj(t) 



(12) 



where Xi is the rescaled GDP defined in eq.© and the 
parameters a(t) and f3(t) can be fixed by imposing that 
the expected number of links 



Lexp(t) = J~)Pt[Xi{t),Xj(t)] 

equals its empirical value [|| 



(13) 



(14) 



and that the expected number of reciprocated links [131 


L Txp(t) = ^Pt[Xi{t),Xj(t)\pt[Xj{t),Xi{t)] (15) 

equals its observed value 0, 13, EH 



(16) 



This particular structure of the WTW topology can be 
tested by comparing various expected topological prop- 
erties with the empirical ones. For instance, we can com- 
pare the empirical and the theoretical plots of vertex de- 
grees (at time t) versus their rescaled GDP (t) 3] . Note 
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FIG. 4: Plot of the rescaled degrees versus the rescaled GDP 
at four different years, and comparison with the expected 
trend. 



that since pt[xi(t), Xj(t)] is symmetric under the exchange 
of i and j, at any given time the expected in-degree and 
the expected out-degree of a vertex i are equal (and this 
is indeed observed in real data, as we now show). We 
simply denote both by k^ xp , which can be expressed as 



p(x) [Til ] . The same is true for any other expected topo- 
logical property [111, [II] • This shows very clearly that the 
WTW topology crucially depends on the GDP distribu- 
tion p(x) shown in fig(21 
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FIG. 5: Cumulative degree distributions of the WTW for 
four different years and comparison with the expected trend. 



k? p (t)=J2pt[xi(t), Xj (t)] (17) 

Since the number of countries N(t) increases in time, 
we define the rescaled degrees fcj(i) = ki{t)/[N(t) — 1] 
that always represent the fraction of vertices which are 
connected to i (the term —1 comes from the fact that 
there are no self-loops in the network, hence the maxi- 
mum degree is always N — 1). In this way, we can easily 
compare the data corresponding to different years and 
network sizes. The results are shown in figsH] for vari- 
ous snapshots of the system. The empirical trends are 
in accordance with the expected ones. Then we can also 
compare the cumulative distribution P> xp \k exp ) of the 
expected degrees with the empirical degree distributions 
P™(k m ) and P™\k out ). The results are shown in figH 
They confirm a good agreement between the theoretical 
prediction and the observed behaviour. 

Note that the accordance with the predicted behaviour 
is extremely important since the expected quantities are 
computed by using only the N GDP values of all coun- 
tries, with no information regarding the ./V 2 trade val- 
ues. On the other hand, the empirical properties of the 
WTW topology are extracted from trade data, with no 
knowledge of the GDP values. The agreement between 
the properties obtained by using these two independent 
sources of information is therefore surprising. Also note 
that all the sums in eqs. (fT3|) . (|T5|) and IT7|) can be rewrit- 
ten in terms of integrals involving only p t [xi (t) , Xj (i)] and 



V. EFFECTS OF THE WTW ON THE GDP 

The above results show that the GDP and its distribu- 
tion determine many topological properties of the WTW. 
As we anticipated, eqs.© and ([7]) suggest that the oppo- 
site is also true, since the GDP is determined by the flow 
matrix representing the WTW as a weighted network. 
Understanding the detailed dependence of the GDP on 
the WTW is an intriguing but difficult problem that is 
currently unsolved. Here we suggest a first step towards 
its solution by proposing a possible framework to be fur- 
ther explored in the future. 

The way the structure of the WTW affects the GDP 
through eq.® is twofold: firstly, its topology deter- 
mines which are the nonzero entries of the adjacency ma- 
trix a,ij(t) and therefore of fij(t); secondly, the nonzero 
weights fij(t) themselves determine Wi(t) quantitatively 
through eq.([3]). From the study of dynamical processes 
on complex networks it is now understood in general 
that the former, purely topological aspect has a primary 
qualitative impact on the collective dynamical behaviour, 
while the latter has mainly a quantitative effect on the 
dynamics. Now we show an example of the crucial role of 
the topology in a paradigmatic case which also turns out 
to be a good candidate for the modeling of our system. 

Interestingly, the effects of a nontrivial network topol- 
ogy on the dynamics of wealth exchange has been ad- 
dressed theoretically in a series of papers [IH, [JJI, [I?], [l8| • 
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In these works, the exchange of wealth is modeled 
through a stochastic process taking place on an under- 
lying network of a gent s. The most general choice which 
has been proposed [17| for such a process is the following 
evolution rule for the wealth Wi of the i-th agent (with 
i = l,N): 

w 2 (t + l) = Wi {t) + &(t) + Vi{t)wi[t) + 

3 

where Wi(t) denotes the wealth of agent i at time t, £i(t) 
is an additive noise term, rji(t) is a multiplicative noise 
term, and Jij(t) represents the fraction of agent i's wealth 
being transferred to agent j. The above model is there- 
fore the combination of an additive and of a multiplica- 
tive stochastic process occurring at discrete time steps 
among the N agents. It is general enough to suspect 
that the combined dynamics of the GDP and the WTW 
can be captured successfully by it. We now discuss how 
the results presented in the literature can be exploited to 
gain insight into this problem. 

For the above model to be fully defined, one needs 
to specify the probability distribution for the stochastic 
variables &(t) and fft(i) as well as the matrix elements 
Jij(t). Several choices have been explored in the litera- 
ture. The most studied case is the purely multiplicative 
one where = and r]i{t) is a Gaussian noise term 
[IM EE Ell- The opposite, purely additive case where 
r)i(t) = and £i(t) is a Gaussian variable has also been 
considered [l7j |. The choices for J-y(t) in all cases assume 
the following dependence on the adjacency matrix (t) 
of the underlying network: 

^) = #)««W = {f (t) !ta$)=0 (19) 

The role of qij{t) is to distribute the wealth coming out 
from i among its neighbours. If all agents deliver the 
same fraction of their wealth to each of their nei ghb ours, 
the choice qij — q/N with constant q is used Il5l [T^ |. 
The choices q^ = q/kj 1 [l6| and q i3 = a I k° ut [PiT| are 
instead considered if the total wealth respectively re- 
ceived or delivered by each agent is a fixed fraction q of 
its wealth, which is equally distributed among its neigh- 
bours. The main point is then the specification of ay(t), 
on which also qij(t) depends. Here the crucial role of 
the topology emerges. Several choices has been explored, 
but in all cases assuming that the topology is held fixed: 
ctij (t) = ciij . In the purely multiplicative case & (t) = it 
was shown that on fully connected graphs (a^ = 1 
Vi,j) the rescaled wealth xi = Wi/{w) eventually ap- 
proaches a stationary distribution displaying a power-law 
tail with exponent determined by q and by the variance 
of the distribution of the multiplicative term rji (t) . Note 
that in this case all the above choices for q^ become 
equivalent since fcj = N — 1 « N. This also implies that 
Jij(t) = q/N The fully connected case is there- 

fore a prototypic example showing that, irrespective of 



the details, the wealth distribution approaches a power- 
law curve. The opposite possibility is that of an empty 
graph where all vertices are disconnected from each other: 
dij = Vi, j. This clearly implies that Jij(t) = Vi, j. In 
the purely multiplicative case, this model is easily shown 
to generate a log-normal wealth distribution since log Wi 
is a sum of independent identically distributed random 
variables, eventually approaching a Gaussian distribution 
as ensured by the Central Limit Theorem. Therefore this 
is the opposite paradigmatic case where, irrespective of 
the model details, the wealth distribution displays a log- 
normal form. 

Interestingly, many empirical wealth and income dis- 
tributions display a mixed shape with a power-law tail in 
the large wealth region and a different behaviour (which 
in some cases is log-normal-like) for the small wealth 
range [l?], EE H3|- This kind of behaviour is also dis- 
played by the GDP distribution across world countries 
as shown in figO Since all real networks fall somewhere 
in between fully connected and empty graphs, it is in- 
teresting to ask whether the observed mixed shape of 
wealth distributions can be accounted for by the topo- 
logical properties of the underlying network. This has 
stimulated the exploration of the model defined in eq. (fl8|) 
in the case of a nontrivial topology [IE, EE EH ■ Remark- 
ably, in the case of purely additive noise on scale-free 
networks [17| and of purely multiplicative noise on het- 
erogeneous networks with varying link density (l8j . the 
wealth distribution has been shown to approach a shape 
very similar to the observed one. 

The above results allow to figure out a likely mecha- 
nism for the GDP evolution driven by the WTW struc- 
ture. Combining eqs. (|H2l3l7p yields the explicit defini- 
tion of GDP in terms of the import/export terms: 

w i (t) = I i (t)+Y / [fji(t)-h(t)} (20) 

3 

It is instructive to compare the above expression with 
the model defined in eq. (fT8|) . First of all, note that the 
GDP is evaluated and publicly released quarterly as the 
main indicator of the overall economic activity of a coun- 
try. Based on the trend exhibited by the GDP at the 
end of each period, each country plans suitable measures 
to adjust its economic activity during the following pe- 
riod. This means that the GDP evolves through discrete 
timesteps in a way similar to the update rule defined in 
eq. (fT8|) . Then, note that the internal demand 7j(t) is an 
endogenous term that does not depend on the GDP of 
other countries, while instead the trade values fij(t) may 
depend on both Wi(t) and Wj(t). Therefore the evolu- 
tion of uii(t) must be the combination of an endogenous 
and an exogenous dynamics. Finally, it is reasonable to 
suspect that, due to the multiplicative character of the 
economy (each country reinvests the income of the pre- 
vious period in the internal activity as well as in the 
trade relationships), the dependece of both the internal 
demand Ii(t) and the trade matrix fij(t) on the GDP 
values is linear as in eq. (fT5|) with rji(t) ^ 0. Therefore we 
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conjecture that a model similar to eq. ([T8|) may capture 
the basic properties of the dynamics of the GDP driven 
by the WTW structure. This possibility will be explored 
in detail in future papers. 

Of course, in a realistic model for the coupled GDP 
and WTW evolution, the strong assumption of a fixed 
topology {flij(t) = ciij) that has been so far used in the 
literature must be relaxed, and the time-dependence of 
the interaction matrix Jij(t) in eq. tfLSf fully exploited. 
On the other hand, we showed that the topology of the 
WTW is at each timestep well predicted by the knowl- 
edge of the values {wi}i as evident from eas. (|12I13I15I17|) . 
Therefore the evolution of (t) is not independent from 
the GDP values, and it should instead contain an explicit 
dependence on them. This dependence should then be 
plugged into the term (t) in eq. (fT5|) . Viewed from the 
point of view of network theory, this property leads us to 
a novel framework where the network is shaped at each 
timestep by the set of fitness values [111 ], which in turn 
evolve in a topology-dependent way through a stochas- 
tic dynamics similar to that in eq. (|18p . The extension of 
the current fitness formalism to take these evolutionary 
aspects into account is an intriguing problem to address 
in the future. 



VI. CONCLUSIONS 

In the present work we have reported a range of 
empirical results and theoretical arguments highlighting 
the interplay between the dynamics of the GDP and 
the topology of the WTW. The topological properties 
of the trade network have been shown to be determined 
by the GDP of all world countries. On the other hand, 



the empirical properties of the GDP distribution across 
world countries can be traced back to the underlying 
dynamical process of wealth exchange among countries, 
which is mediated by the WTW. The emerging picture 
is that of a discrete process where at each timestep 
the GDP distribution is determined by the WTW as 
a weighted network, and at the same time the WTW 
topology is also determined by the GDP values. We 
have thus suggested the need for a theoretical framework 
where the network topology is determined by some hid- 
den quantity, which is in turn an evolving variable (and 
not simply a quenched one) whose dynamics is coupled 
to the topology. The present work provides robust 
empirical evidence for such a framework, highlighting 
the need to further develop current models in order 
to take evolutionary aspects into account. We have 
proposed a paradigmatic class of stochastic models for 
the fitness variable that may suggest a possible way to 
integrate the dynamics with the topology. The present 
work represents a basis for a future understanding of 
the details of this interplay in the case of the world 
economy. Such an improved framework would give 
predictive results on extremely important issues such 
as the economic power and interdependence of world 
countries, but also on many other problems in social 
science and biology where formally similar mechanisms 
are at work. 
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